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A  FREE  ELECTRON  LASER  WITH  A  ROTATING  QDADRUPOLE  WIGGLER 


I.  Introduction 

1  2 

Free  electron  lasers  *  are  devices  which  are  designed  to  produce 
coherent  radiation  by  passing  a  relativistic  electron  beaa  through  some  sort 
of  spatially  periodic  perturbation.  Commonly,  a  periodic  magnetic  field, 
called  a  wiggle r,  is  used. 

In  general,  the  wlggler  field  varies  transversely  as  well  as  axially. 
However,  in  most  experiments  the  initial  beam  radius,  R^  is  a  small  fraction 
of  the  wlggler  period,  \ and  the  transverse  variations  of  the  wlggler  are 
neglected.  For  example,  the  field  of  a  periodic  helical  winding  of  two  wires 
around  a  drift  tube  to  lowest  order  in  c  ■  k^R^  (k^  ■  2w/\^)  Is  commonly 
represented  by  a  simple  periodic  function  of  axial  distance  alone 

"  B^lcos  +  -lB  (!) 

We  provided  the  notation  of  this  field  with  subscript  "d"  in  order  to  indicate 
that  such  a  field  has  a  dipole  character. 

Generally,  a  second  field  in  the  axial  direction  is  added  to  confine  the 
beam  of  particles  against  their  natural  repulsion. 

In  principle,  one  can  generate  a  periodic  magnetic  field  by  a  helical 

winding  of  1%  wires  (1  ■  1,2...).  A  magnetic  field  of  this  kind  exhibits 

helical  symmetry  and  it  is  convenient  to  describe  it  by  means  of  a  scalar 

potential  that  satisfied  Laplace's  equation.*  This  potential  depends  only 

on  dimensionless  coordinates  p  ■  a  r  and  u  ■  fl  -  a  z  and  has  a  fora 

w  WWW 

♦  ,  ■  — -  I£(tp)  tu  •  (2) 

_  w 
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Where  1^  is  a  modified  Bessel  function  of  order  £,  the  pitch  of  the  helix 
is  defined  by, 

<*w  -  v*  (3) 

and  8  is  the  polar  angle.  By  definition  the  components  of  the  corresponding 
magnetic  field  can  be  deduced  from 

£  -  grad  .(4) 

Consider  now  that  the  periodic  magnetic  field  is  generated  by  a  helical 
winding  of  four  wires  as  shown  in  Fig.  1  (the  current  flow  is  the  same  as  in  a 
stellerator  configuration) . 

We  expand  l£  to  second  order  in  e  and  from  Eq.  (4)  we  obtain 

Br  a  ■  Bno  k«r  Si°  <29  ~  k  *> 

r,q  oq  q  q 

V,  '  ®oq  V  «•  <“  ‘ 

■  hW'£)2  co*  <2e '  V>  <5> 

where  BQq  is  the  value  of  the  magnetic  field  at  r  -  X^/2w ,  and  the  subscript 
q  indicates  that  the  field  is  quadrupole  in  nature.  To  lowest  order  in  e  from 
Eqs.  (5)  we  obtain  for  the  magnetic  field  components  in  cartesian  coordinates 


B  k  (y  cos  k  z  -  x  sin  k  z) 

oq  i  q  q 


k  (  *  cos  k  z  +  y  sin  k  z) 
oq  q  q  q 


B  -  0. 

zq 


006  can  show  chat  Che  field  lines  of  Che  above  magnetic  field,  Eqs.  (6),  are 
similar  to  Che  line  forces  of  a  quadrupole  magnet  continuously  rotating  along 
the  z  axis.  Quadrupole  fields  are  known  in  accelerator  physics  for  their 
focusing  ability  of  the  beam  current.^ 

In  this  paper,  we  investigate  the  feasibility  of  using  a  quadrupole 
magnetic  field  as  a  wlggler  in  free  elelctron  laser  devices  as  well  as  a 
focusing  field  to  confine  the  beam.  Ue  assume,  that  if  «  1,  the  actual 

field  can  be  good  represented  by  Eqs.  (6).  Because  of  the  strong  focusing 
property  of  the  field,  one  might  expect  that  large  beam  currents  could  stably 
propagate.  However,  because  the  field  is  zero  on  the  axis,  the  free  electron 
laser  mechanism  might  be  weaker.  Here  we  investigate  these  Issues. 


II.  Particle  Orbits  and  Stability  Analyses 

The  physical  model  we  develop  will  consider  a  relativistic  non  neutral 
electron  beam  with  radius  Rb  propagating  in  the  quadrupole  magnetic  field 
given  by  Eqs.  (6).  We  assume  that  the  beam  is  solid  and  has  a  uniform  density 
nb(r)  -  nQ.  The  expressions  for  the  equilibrium  self-electric  and  self- 
magnetic  fields  associated  with  the  beam  space  charge  and  axial  current 
Jb  ■  “«n0vb  (vjj”beam  axial  velocity)  are  given  by. 


E  -  -  2ir|e|n  x 
xs  o 

E  «  -  2»|e|n  y 
ys  oJ 

®X8  “  2*|e|Vby> 

and 

Bys  “  "  2^lelVbX*  (7) 

for  0  <  r  -  (x2  +  y2)1/2  <  1^ 

where  |e|  Is  Che  electron  charge  and  0,  ■  v./c.  The  electron  orbits  within 

D  D 

the  beam  are  determined  from  the  equations  of  motion 

df<YVx)  "  I  “ b(1  "  W*  +  Vxnqkq(x  C08  kqZ  +  7  Sin  kq° 
df  iTfy)  "  7  “b(1  '  0*6b)y  +  v*°qkq<X  8in  kqZ  "  7  C08  kqz) 

df  (YVz}  "  I  “b8b(x8x  +  78y>  +  flqV  V7  C08  kyz  "  x  8ln  kqz) 

-  v  (  x  cos  k  z  +  y  sin  k  z)l  (8) 

x  q  q 

2  2 

where  co^  -  AimQe  /m  Is  the  beam  plasma  frequency  and  •  lelB^/mc  is  the 
cyclotron  frequency.  From  the  above  equations  and  the  assumption  of  this 
analysis  (k  x  (k  y)  -  e  «1)  is  evident  that  v  /v  *>■  c,  and  v  /v  -  e, 


therefore  to  lowest  order  in  e  the  system  (8)  reduces  to  a  simple  system  of 
equations  describing  the  transverse  motion  of  an  electron. 


2 

dvx  j  <n£  n 

Vz  17  “  2  T  *  +  v*  71  V*  CO*  V  +  7  8in  V> 

Yo  0 

2 

dv  i  a.  a 

Vz"df"2“37  +  VZ^  kq(x  Sln  kqZ  “  y  co8  kqz)  («> 

7o  0 

where  T  -  Y  *  const,  and  8.-8  a  const.,  hence,  ~  -  v  ~. 

o  b  z  dt  z  dz 

We  Introduce  dimensionless  variables;  the  distance  is  multiplied  by  kq 
and  frequency  is  divided  by  (kqc).  Then  the  normalized  Eqs.  (9)  have  a  form 


— f  *  a  (*  cos  z  +  y  sin  z)  +  Sx 
dzz 


■  a  (*  sin  z  -  y  cos  z)  ♦  Sy 


(10) 


where 


a  h  nq/(yo0zkqc) , 

6  5  *S'<2yokqc2>  (11> 

a  measures  the  strength  of  the  quadrupole  fields  and  8  measures  the  beam 
density.  Although  we  did  not  change  the  notation  for  x,  y  and  z,  one  should 
regard  these  variables  as  dimensionless. 

In  order  to  solve  the  system  of  Eqs.  (10)  it  is  convenient  to  introduce 
new  variables  u  and  v  related  to  x  and  y  by  the  following  transformation 


x  ■  u  cos  (*/ 2)  +  v  sin  (z/ 2) 


y  ■  u  sin  (z/2)  -  v  cos  (z/2). 


(12) 


In  the  new  variables,  Eqs.  (10)  becomes  a  system  of  homogeneous  equations  with 
constant  coefficients,  namely 


d2u  +  dv 
dz2 


dz 


dfv 

dz2 


du 

dz 


(a  +  I  +  6)  u 
4 

(a  -  7  6)  v 

4 


0 


0. 


(13) 


The  solution  of  this  system  of  equations  has,  in  general,  the  following  form 


u  »  A  exp  ikz 


and 


v  ■  B  exp  ikz.  (14) 

Inserting  this  solution  into  Eqs.  (13),  we  obtain  a  system  of  two 
algebraic  equations  to  determine  the  coefficients  A  and  B.  Non-triviality  of 
the  solution  of  that  algebraic  system  provides  the  equation  for  k 

(k2+  ^  +  5)2  -  a2  -  k2  -  0.  (15) 

The  four  Independent  solutions  of  the  bi-quadratic  equation,  (15),  are  given 
by 


6 


If  one  of  the  k8  has  a  non-zero  Imaginary  part  the  solutions  (14)  become 
unstable.  Therefore,  for  I*  kg  -0,  the  following  conditions  are  required 

/J  <  a  <  («  +  1/4)  (17a) 


or  equivalently. 


2  2Toaa 


,,  2  .  .  2  2  3.2 

(2u.  +  k  c  y  ) 

D  q  o 


e 


3  2  2 

c 

o  q 


(17b) 


In  the  limit  of  negligibly  small  self  fields  (5  ♦  0)  the  stability  conditions 
(17)  become 


a  <  7* 


Equation  (17)  indicates  that  a  quadrupole  magnetic  field  should  be  strong 
enough  to  overcome  the  repulsive  forces  due  to  self  fields  of  the  beam,  but 
not  too  strong. 


The  stable  region  in  the  parameter  space  or  and  6  is  shown  in  Fig.  2 

2  2  2 

between  the  two  solid  lines,  a  ■  6  and  a  *  (6  +  1/4)  .  It  is  convenient 
relate  the  parameter  6  to  the  beam  current*  The  relation  is 


V1.  ■ ,n  (#.v2  <k,v2 


where  I  *  me  /e8_Y-  is  the  Alfven-Lawson  current. 


From  Eq.  (17)  it  follows  that  the  theoretical  maximum  for  6  In  the  stable 

2 


region  is  determined  by  5  *  (d  +  1/4)  and  is  equal  to  1/4.  Therefore,  In 
principle,  the  maximum  electron  beam  current  which  is  able  to  propagate  in  a 
quadrupole  magnetic  field  is 


wK  <wJ 


I  . 
a 


(18) 


For  example,  if  Y  ■  4,  k  R.  *  7,  1  “3  cm,  I,  a  8.5  kA,  which 

o  q  b  a  q  b,max 


require  B  *  3.5  kG. 
oq 


It  is  also  interesting  to  compare  the  maximum  beam  density  which,  can  be 
confined  by  a  rotating  quadrupole  field  with  that  which  can  be  confined  by  a 
uniform  axial  magnetic  field.  It  is  not  difficult  to  show  in  the  latter  case 
that  the  condition  for  stable  orbits  is 


2  YoQc 


(19) 


where  Q  is  the  cyclotron  frequency  in  the  axial  field.  Comparing  Eq.  (19) 


with  Eq.  (17b),  it  is  apparent  that  the  quadrupole  field  with  -  ftc  can 


confine  four  times  the  beam  current. 


III.  Single  Particle  Analyses  of  a  Free  Electron  Laser  with  Quadrupole  Wlggler 

In  the  previous  section  we  investigated  the  stability  of  the  electron 

orbits  in  a  relativistic  electron  beam  propagating  in  a  quadrupole  magnetic 

field,  Eqs.  (6).  Here,  we  will  show  that  such  a  field  can  be  used  as  a 

wlggler  in  a  free  electron  laser  device  to  produce  coherent  radiation  with 

2 

wavelength  A  *  A  /2v  .  Henceforth,  a  subscript  r  denotes  the  radiation 


In  our  analyses  we  will  use  a  single  particle  approach  to  describe  the 
interaction  of  a  relativistic  electron  with  the  periodic  quadrupole  magnetic 
field  and  a  plane  circulary  polarized  electromagnetic  wave.  This  is  analogous 
to  the  calculation  of  Colson^  of  free  electron  laser  gain  in  a  conventional 
wlggler  (Eq.  (1)),  oscillator  configuration  with  a  specified  radiation 
field.  The  radiation  is  represented  by  a  circularly  polarized  electromagnetic 


E  [cos  (k  z  -  to  t  +  <b)i  -  sin  (kz  —  u>  t  +  <fr)  i  I 
ol  r  r  y  x  r  r  y  yJ 


B  -  i  x  E 
— r  z  — 1 


with  ix,  iy  and  i2  being  unit  vectors  in  x,  y  and  z  direction,  respectively, 

and  (d  ■  k  c. 
r  r 

The  equations  of  motion  for  the  relativistic  electron  in  the  combined 
electromagnetic  wave  (20)  and  pump  wave  (wiggler)  (6)  are  given  by 


d  dx 
dx  Y  dx 


-e  (1  -  0  )  +  0  b 
rx  z  z  qy 


■p  Y  -  e  (1-0  )  -  B  b 

dx  dx  ry  z  z  qx 


__d  dz 
dx  ^  dx 


-  (0  (b  -  p  )  -  0  (b  +  p  ). 
y  qx  ry  x  qy  rx 


(21a) 


(21b) 


(21c) 


As  before,  we  kept  only  terms  of  lowest  order  in  e  and  we  introduced 
normalized  variables.  Ve  also  denoted 


■  era[cos  (x(z  -  x)  +  $}ix  -  sin  (k(z  -  x)  +  $!iy  (22) 


9 


where 


b  *  b  !(▼  cos  z  -  x  sin  zli  +  fx  cos  z  +  y  sin  z!i  I  (23) 
-q  qa 1  J  ’  x  yJ 


et.  *  W  V<*CV) 

-  v'V’  «*) 

and  tc  ■  k  /k  ,  t  -  (k  c)t  x,  y  and  z  are  also  dimensionless  (note  that 
r  q  q 

ic  «  2y2  »  1). 
o 

The  orbit  equations  (21)  describing  the  motion  transverse  and  parallel  to 
the  electron  stream  are  coupled,  due  to  the  x  and  y  dependence  of  the  pump 
field.  This  fact  makes  direct  integration  difficult,  even  making  a  linear 
approximation.  The  principal  complication  is  that  the  orbit  has  two  very 
different  frequencies,  a  forced  oscillation  at  the  doppler  shifted  pump  or 
radiation  frequency,  and  a  slower  evolution  at  the  beat  frequency.  To  proceed 
we  will  average  over  the  fast  oscillation  and  calculate  the  motion  only  on  the 
slow  time  scale. 


IV.  Multiple-Time  Scale  Analyses 


In  order  to  find  the  solution  of  the  set  of  the  nonlinear  coupled 


equation  (21),  we  will  utilize  multiple-time  scale  analyses.  This  gives  the 


same  result  as  a  similar  derivation  making  use  of  the  covariant  structure  of 

g 

the  equations  of  motion.  We  introduce  two  time  scales,  one  fast  -  t  , 

o 

another  slow  -  Tj.  The  fast  time  scale  is  associated  with  the  pump  field 
"frequency"  -  (k^c)  1 ,  and  the  slow  time  scale  is  proportional  to  k  v^ , 
which  in  turn  is  proportional  to  0^  1  and  (|e|Eo/rac)  l.  Also,  in  the  slow 


10 


■*.  % 


time  equations  we  retain  as  a  low  frequency  term  the  beat  frequency 
\  -  Ckr  *  VV 

We  assume  that  variations  on  fast  time  scale  are  periodic,  and  both 
variables  to  and  are  treated  as  independent  variables.  To  introduce  the 
difference  between  fast  and  slow  time  scales  we  consider  a  formal  procedure 
consisting  of  assuming  a  perturbation  expansion  of  the  form 

x( r© , x i )  »  x0^T0»Tp  +  VVV  +  *2*To,TP  +  *** 


and 


_d  _3_  _3_ 

dt  "  3t  3t, 
o  1 


(25) 


from  which  follows 


32 

S'?7! 


where  we  denoted  the  order  of  the  quantities  with  a  subscript  o, 1,...,  hence 
3  3  32 

*o  an<*  xi  TZ-  an<l  x2  an<*  - J  are  *®fo,  first,  and  second  order 

o  aTl  3t. 

quantities,  respectively. 

Substituting  Eqs.  (25)  into  Eqs.  (21)  for  zero  order  equations  we  obtain 


3 


3x 


o  3t 


0 


3  3yo 

TT\TT 

o  o 


0 


3z 


'o  3t 


0. 


(26) 


11 


The  solutions  to  Eqs.  (26)  are 


VvV 


WV  - 


xo(Tl) 

*o<V 


•o^l’V  "  zo<Tl>  +  8o(Tl)To  (27) 

where  3z  /3t  ■  8  (t.)  and  3x  /3t  ■  3y  /3t  -  0,  since  we  assumed  that  zero 

oo  oi  oo  oo 

order  transverse  motion  doesn't  exist.  The  quantity  as  expressed  in 

Eq.  (27)  is  undefined  because  any  velocity  which  is  nearly  equal  to  the 

particle  velocity  will  render  z^ij)  «nd  slowly  varying  functions  of 

time.  To  be  more  specific,  we  define  8  (t.)  -  v  /c  in  terms  of  the  resonant 

oi  o 

velocity,  so  that 


v 

o 


k 


r 


(28) 


Since  no  quantity  on  the  right-hand  side  of  Eq.  (28)  depends  on  time,  80(Tj) 
is  independent  of  tj.  He  will  now  show  that  this  is  consistent. 

The  first  order  part  for  Eq.  (21c)  is  given  by 


if  <*.  if>  ♦  if  (V„>  *  if  (Vo>  *  if  <T.  if>  ■ 

O  O  1  0  O  1 


(29) 


Upon  averaging  over  the  fast  time  scale  we  have. 


(r,B0) 

o  o 


0. 


(30) 
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T 


T 


Note  that  y  depends  only  on  8  ,  therefore  Eq.  (30)  requires 
o  o 


8  (t,)  ■  const. 

O  i 


Thus  our  assumption  that  8  is  Independent  of  t,  is  consistent.  Making  use  of 

o  1 

this  fact  we  return  to  Eq.  (29)  and  Integrate  it  once  in  tq.  The  resulting 
equation  is  given  by 

hr- *  *  V. '  c<Ti>  <31> 

o  1 

where  C(t^)  is  a  constant  of  the  Integration. 

Expanding  y  to  the  first  order  and  utilizing  Eq.  (27) ,  we  find  that 


V 


1 


y3b 
o  o 


(32) 


Substituting  it  into  Eq.  (31)  we  obtain 


C(t.) 


■  0. 


(33) 


The  last  two  terms  do  not  depend  on  tq,  therefore  by  averaging  Eq.  (33)  over 

t  we  obtain 
o 


Substituting  this  result  back  into  Eq.  (33)  we  find  that 


o 


13 


VS 


follows,  chat 


W  *  Vo8..!'1!* 


where  3z  /3t  «  0  ,(t.). 

O  X  Z ,  1  X 

Consider  now  the  first  order  transverse  motion.  The  left-hand  side  of 
the  Eq.  (21a)  becomes 


.  3x  „  3x  3x  3x, 

3  y  _2.  +  JL  f  Y  _ 2.  +  v  _ 2.1  +  _£_v  __1 

3t  'o  3t  3t  ^To  3t.  3t  '  3t  ro  3r 

o  o  O  1  0  0  o 


Since  xQ  does  not  depend  on  To>  the  above  expression  reduces  to 


3t  To  3t  * 
o  o 


(35a) 


and  correspondingly,  for  left-hand  side  of  Eq.  (21b) 


3t  yo  3t 
o  o 


(35b) 


The  radiation  and  wlggler  field  expressions  have  the  following  fora 


*r,I  ‘  *r,[co*  >'!<8o  ‘  1,To  +  'o'  +  ♦>*, 


-  sin  (<(S0  -  I)  to  +  so)  +  +|iy,] 


V  .V 


The  dependence  of  the  right-hand  side  on  Is  through  zQ  which  Is  a  function 
of  tj.  We  proceed  with  the  multi-time  expansion  procedure.  To  do  so.  we 
write  the  second  order  equations  for  Eqs.  (21).  The  expressions  for  the 
wiggler  field,  evaluated  along  the  first  order  orbit,  and  therefore  correct  to 
second  order  are  given  by 


b  -  b  l  y  cos  (z  +  It)  +  y,  cos  (z  +  B  t  )  -  z  y  sin  (z  +  B  t) 
qx  qa 1  o  o  oo  1  o  oo  1  o  o  oo 

-  x„  sin  (2  +  B  t)  -  X,  sin  (z  +  It)  -  z.x  cos  (z  +  It  )]  (42) 

O  OOO  1  OOO  10  o  o  o  1 


b  «  bl x  COS  (z„  +  B  t  )  +  x.  cos  (z  +  It)  -  z.x  sin  (z  +  8  T  ) 
qy  qal  o  o  00  l  o  00  io  o  00 


+  yo  sln  <zo  +  Vo*  +  yl  sln  (zo  +  Vo*  +  zlyo  co8  (zo  +  BoTo)l*  <43* 


However,  we  are  interested  in  deriving  equations  governing  the  particle  notion 
on  the  slow  tine  scale  t^.  Therefore  we  perform  an  average  over  the  fast  tine 
scale.  We  denote  this  averaging  process  by  <  >  brackets,  for  example,  Eq. 
(20a)  becomes 


<(•57  7  Hv  •  -  <1*0  «  -  W  ♦  «»zV>2> 


qy 


(44) 


2  2 

For  the  left-hand  side  of  Eq.  (44)  after  averaging  we  obtain  (3  xq/3Tj  ), 
Thus  the  second  order  equation  becomes 


32x 


■  <e*,l(*r*,l  +  b,y.l»  '  <*rx,2>  +  So  <en>,2  *  V,2>- 


a  t 


(45) 


Inserting  Eq.  (43)  into  Eq.  (45),  using  the  fact,  that  ^  and  do  not 
depend  on  tq  and  the  assuaption,  that  3erx/3x(}  ■  Se^/dy^^  -  0,  we  obtain 


r0  7T  "  8obq.  <XI  cos  (*o  +  aoTo)  +  yI  sin  (zo  +  8oTo)>-  <46> 

3T1 


Analogous  to  Eq.  (46)  we  can  write  an  equation  for  the  y~conponent 


T  - 7  “  ®  b  <*,  (*  ♦  8  T  )  -  y.  cos  (z  +  B  r  )>.  (47) 

o  3  2  oqal  ooo  1  ooo 

Making  use  of  the  expressions  for  x^  and  y^  (Eqs.  (40)-(41))  end  averaging 
over  the  fast  tlae  scale  we  obtain 


3*x  e  b 

- r  ■  ■  cos  | (<  +  I)x  ♦  4j 

3tI  V  ° 


-Jix 

'l  ° 


3  y  e  b 
Jo  _  ra  qs 


!<«♦  »>«o  *  *J 


In  writing  Eqs.  (48),  we  have  utilized  the  resonance  condition,  Eq.  (28), 

which  here  is  expressed  as  8  ■  x/(x  +1).  Consider  now  the  second  order 

o 

equation  for  Eq.  (21c).  The  left~hand  side  of  Eq.  (21c)  after  averaging 
3 

becomes  Y0(3$zj/3Tj),  here  we  also  used  relation  (34).  Thus  we  can  write. 


J  j2'o 


7T?  ■  <e7.><h«y.i  •  V.i>  '  *  V,i»-  <*» 


To  obtain  the  expressions  for  0  ,  and  0  ,  we  differentiate  Eqs.  (40)-(41), 

7  9*  Xg 1 

respectively.  Then  utilizing  Eqs.  (36)— ( 37 )  and  averaging  over  the  fast  time 


scale  we  obtain  the  z-component  of  the  equation  of  action  on  the  slow  tlae 
scale,  naaely 


.  ILM  J 


Y  8 
To  o 


COS  l  (if  +  l)zQ  +  *] 


-  XQ  sin  1<K  +  1)zq  +  ♦  ]}. 


(50) 


(Recall  that  zQ  is  not  the  unperturbed  orbit,  but  is  the  orbit  minus  the 
resonant  velocity  tines  tine,  see  Eq.  (27).) 


Denoting 


a  _  2  * 


Z 

Yo*c 


Z  Y^S 
'opo 


a  = 

a  -  v  • 

To 


(51) 


the  final  set  of  equations  of  motion  on  the  slow  time  scale  Is 


d2x 

2  ” 

a 

cos  [ (k 

+  1)* 

+  $J  -  a2x 

(52a) 

dx 

d2y 
— 2  “ 

a 

sin  l (k 

+  1)2 

+  ♦]  -  ®2y 

(52b) 

dt 

y  cos 

K 

K  +  1)* 

+  ♦]  * 

x  sin  [  (pc  +  l)z  +  *  J } . 

(52c) 

For  convenience,  we  have  dropped  the  subscripts  o  and  1,  since  for  the 
following  analyses  we  will  use  only  the  Eqs.  (52).  When  necessary  we  will 
return  to  the  original  notation. 


V.  Linear  Analyses 

In  order  to  solve  the  non-linear  coupled  system  of  Eqs.  (52)  we  will 
employ  a  perturbation  analysis.  We  assume  that  EQ  (a,  az)  la  small  and  the 
electron  displacement  and  energy  change  can  be  expanded  in  powers  of  EQ. 

For  convenience  of  the  analysis  we  introduce  a  complex  variable 
£  ■  x  +  iy,  in  which  case  Eqs.  (52)  become. 


a  exp  (ix)  -  a  £ 


.2  a 

— 2  "  21  U  exp  (“  ix)  -  K  exp  (ix)] 
dr 


(53) 

(54) 


%  W 

where  x  *  (*  +  l)z  +  ♦  =  <z  +  $  and  £  is  the  complex  conjugate  of  C*  The 
equations  to  zero  order  in  the  radiation  field  are 


dVo) 

dT2 


_  2  (o) 

«  t  » 


dV»> 

^7“ 


•  0. 


(55) 


The  general  solution  of  the  first,  harmonic  oscillator,  equation  is  given  by 


r(°)  _  i(<*T  +  0  )  -i(crr  +  0  ) 

;  P+e  +  +  p_e 


(56) 


where  the  real  quantities  p^,  p_,  0+  and  0_  are  obtained  from  the  Initial 


conditions.  In  Appendix  A  we  derive  these  quantities  in  terms  of  the 
parameters  describing  the  injection  of  the  electron  beam  into  the  wiggler 
field.  The  oscillation  corresponds  to  betatron  oscillations  in  the  strong 
focusing  wiggler  field. 

The  solution  of  the  second  equation  in  (55)  is 

-  *1  +  (8£  -  6q)t  (57) 

where  z^  and  8^  are  the  position  of  the  electron  along  the  z  axis  and  its 
velocity  at  the  moment  (t  »  0)  of  entering  the  interaction  region,  - 
respectively.  (Recall  that  z^°^  in  Eq.  (57)  varies  slowly  in  t.  The  actual 
fast  variation  is  obtained  by  adding  a  (®0T)  to  z^°\) 

Ue  denote 

Aw  •  k(8.  -  8  ) 

1  o 

using  the  resonance  condition  we  can  reexpress 

Aw  -  o>r  -  (kr  +  kq)v£. 

Thus  the  zero  order  expression  for  x  becomes 

X<0)  -  *0  +  Awt  (58) 

where  ♦  ■«.+*. 

O  1 

We  now  calculate  C  to  first  order  in  E  .  It  obeys  a  forced  oscillator 

o 

equation  with  ■  d£^/dt  »  0  at  t  -  0  as  initial  conditions,  so  it  has  a 


a  exp  (i*o)/(a+a_){  exp  (iAwx)  -  —  exp  (iax) 


"  I  T  exp  (“  iaT)f 


where  *  a  ±  Aw. 

The  first  order  expression  for  z  is  obtained  by  integrating  Eq.  (54)  with 
C  ■  5^°^  and  x  *  X(0)  »  th*  result  is  given  by 


z  *  a 


,{ — j  [sin  (a+x  +  *Q  +  0_)  -  sin  (0_  +  *q)  -  a+r  cos  (0_  +  *Q)j 


- r  [sin  (a  x  —  ♦  +  0.)  -  sin  (0.  -  *  )  -  a  x  cos  (0,  -  *  )]}.  (60) 

4  w  “  O'  t  o  +  o 


The  equation  for  the  second  order  term  in  the  expansion  of  z 


d2z(2)  az  r  (l)|r(o)  f  .  (o).  .  r(o)*  fI  (o). . 

- 2 —  m  -  —  \<z  [C  exp  (-  iy  )  +  5  exp  (ix  )j 

dx 


Mc<1).xp(-ix(o,)-5(I,*.xp(ix<'»)|t. 


The  first  two  terms  on  the  right-hand  side  are  similar  to  those  of  a 
conventional  free  electron  laser,  in  that  only  zero  order  transverse  motion  in 
the  wlggler  comes  into  play.  However,  the  frequency  of  this  term  is  shifted 
due  to  the  betatron  oscillations  of  the  electrons.  The  second  two  terms  on 
the  right-hand  side  have  no  analog  in  the  standard  one-dimensional  treatments. 
They  describe  a  resonant  beating  of  the  first  order  transverse  oscillation 
with  longitudinal  motion.  These  last  two  terras  on  the  right-hand  side  of  Eq. 


(61)  afCer  Inserting  Eqs.  (58)-(59)  will  take  the  following  form 


a  a 
2 

-= -  (a  sin  a.T  -  a.  sin  a  r). 

2aa+a_  -  +  + 


Inserting  Eqs.  (58)  and  (60)  into  the  first  two  terms  in  the  right-hand  side 

of  Eqs.  (61)  will  result  in  an  expression  which  will  include  terms  dependent 

on  <pQ.  We  assume  that  particles  are  injected  uniformly  in  z,  so  that  an 
average  over  can  be  done.  Thus,  the  average  expression  for  the  right-hand 
side  is  given  by 

2  - 

a-  *  P+P_  1  1 

- 2 -  J - jp)  8itl  ^2aT  +  ®+  +  9  ^ 

a.  a+ 

-  Ip+  COS  (o_T  +  6+)  +  p_  COS  (o+r  +  0_)J 

p+  p- 
X  [— J  (sin  0+  +  o_T  cos  0+) - 2  +  a+T  cos  9_)] 

a-  ’  a+ 

+  [p+  sin  (o_x  +  0+)  -  p_  sin  (a+r  +  0_)] 


p+  p- 
x  [— j  (cos  0+  -  a_r  sin  0+)  +  — ^  (cos  0  -  a+x  sin  0  )]. 

3-  " 

At  that  point  we  introduce  an  additional  assumption  by  taking  p_  »  0, 
this  simplifies  the  last  expression  considerably.  In  Appendix  A  we  show  that 
this  assumption  represents  initial  conditions  for  an  electron  whose 
cylindrical  radius  remains  constant  in  time  and  rotates  in  azimuthal  angle 
with  angular  velocity  a.  An  analysis  utilizing  more  general  initial 
conditions  will  be  performed  in  the  future. 

Thus  Eq.  (61)  reduces  to  the  following 
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2  (2)  a2<pl 

G  Z  Z  '  /  ,  v 

- 5 —  ■  - 5—  (sin  a  t  -  a  t  cos  a  t) 

d-r  2a  1 


a  a  sin  a.T  sin  a  t 

+  -A-  ( - ± - —  1. 

2a  v  a+  a_ 


Integrating  Eq.  (62)  once  we  obtain 


/-»  a  a  1  -  cos  a.T  1  -  cos  a  t 
2a  [  2  2  ' 


azKp+ 


(1  -  cos  a_T  -  y  a_i  sin  a_i). 


VI.  Single  Pass  Gain  Calculation 

Following  the  usual  procedure**  the  gain  G(t)  is  defined  by 


Y.  -  y(t)  ,  n  V 

G(t)  -  - ymc^-f-  (64) 

Yi  1  Ez  V 

o 

47 

where  V  is  the  volume  of  a  section  of  the  beam  with  number  density  ng, 

_  2 

Y(t)mc  is  the  average  over  of  the  electron  energy. 

The  averaging  of  y  *  y  +  yf^  +  yf^  +  •••  over  #  eliminates  the  first 

Oil  o 

order  contribution  in  (64),  therefore  inserting  Eq.  (63)  into  y,  ^  ■  y^p  6'^} 

*  1  o  o  z ,  1 

and  the  result  into  Eq.  (64)  yields  the  following  expression  for  the  gain 


a_a  (1-cosa  t)  (1-cosa.r) 


-3 —  (1  -  cosa_x  -  ~2  (ajr)sin  a_t)f, 


Utilizing  the  definitions  for  a,  a  ,  a  and  6  we  obtain 


(1-cosa  t)  (1-coso  t) 


G<T)  "  i#«  1  (a  t)2 


( a+T  ) 


-1} 


-  {p2Y2  (1  -  cosa_T  -  y  a_T  sina_x)/(a_x)3f 


(66) 


where 


G(x)  -  G(x)/t4(-f)  a2(ax)3]. 


We  denote  the  first  curly  bracket  in  Eq.  (66)  G^(r)  and  the  second  Gp(x),  thus 

AAA  * 

G  ■  +  Gp.  The  expression  for  Gp(x)  is  similar  to  the  formula  for  gain 

obtained  in  single  particle  analyses  with  a  dipole  wiggler  (see  Eq.  (16)  in 
Ref.  (6)).  The  main  difference  is  that  the  resonance  is  at  Aw  *  o  instead 
of  Aw  -  0.  Thus  the  beat  frequency  resonates  with  the  betatron  frequency. 

Had  we  chosen  instead  the  other  initialization,  p_  4  0,  p+  ■  0,  p_  *  0,  the 
particles  would  rotate  in  the  opposite  sense  and  the  resonant  condition  would 

A 

have  been  Aw  *  -  a.  Gfe  is  an  additional  term  associated  with  the  transverse 
bunching  in  the  FEL  with  a  quadrupole  wiggler,  however  it  is  smaller  in 

A  A  2 

magnitude  than  Gp(x).  In  Fig.  3  we  show  the  (G^/S*)  as  a  function  of  (ax) 

*  2  2 

for  different  ax  and  in  Fig.  4  we  present  (G  /p^Y  )  as  a  function  of  (a  x). 

p  +■  o  - 

A 

A  plot  of  the  normalized  gain  is  shown  in  Fig.  5.  The  maximum  of  Gp(x) 
appears  at  a_x  *  2.6  and  is  given  by, 

V„x  *  °'27 


This  is  nearly  equal  to  the  maximum  value  of  G.  Defining  t  «  L/c6q,  where  L 
is  the  interaction  region  length,  the  maximum  gain  is  given  by 


G 

max 


0.27  (e4B2  n  X  )( 
oq  e  q 


2*R, 


V 


(67) 


where  N  -  L/\  is  the  number  of  wiggler  periods  in  the  interaction  region. 

<1 

From  Eq.  (67)  we  see  that  G„av,  for  an  FEL  with  a  quadrupole  pump  has  the 

*  UiaX 

same  functional  dependence  as  an  FEL  with  a  dipole  pump,  provided  ®oq^q^b^‘ 


Conclusions 

We  have  shown  that  relativistic  electron  trajectories  in  a  rotating 
quadrupole  periodic  magnetic  field  are  stable  even  if  large  space  charge 
forces  are  present.  This  is  not  the  case  for  the  electron  trajectories  in  a 
dipole  periodic  magnetic  field.  In  the  last  case  to  improve  the  stability  a 
constant  magnetic  field  is  introduced  in  addition  to  the  periodic  field. 
However,  analyses  have  shown  that  the  orbits  are  not  stable  for  all  Initial 
conditions.^* In  the  case  of  a  quadrupole  field  the  stability  conditions 
and  the  theoretical  maximum  value  for  the  electron  beam  current  depends  on  the 
beam  energy,  density  and  radius  and,  in  principle,  can  be  close  to  the  Alfven- 
Lawson  limiting  current. 

The  linear,  low  gain  analysis  indicated  that  by  utilizing  a  quadrupole 
periodic  magnetic  field  as  a  wiggler  in  an  FEL  we  obtain  a  positive  gain.  The 
maximum  value  of  this  gain  has  the  same  parametric  dependence  as  for  an  FEL 
operating  with  the  usual  dipole  wiggler,  provided  the  same  value  of  the 
amplitude  of  magnetic  field  at  the  beam  position  can  be  produced. 

Therefore,  an  FEL  with  rotating  quadrupole  pump  represents  an  interesting 
new  concept  to  obtain  high-power,  coherent  radiation  in  the  millimeter  and 
sub-millimeter  regime. 

In  a  future  work  we  will  investigate  the  feasibility  of  such  devices  by 
performing  single-particle  nonlinear  analyses  in  Compton  regime  and  an 
analysis  of  the  device  in  Raman  regime. 
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directions  to  produce  a  continuously  rotating  quadrupole  magnetic 
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Shows  the  stability  region,  where  the  upper  line  is  a  -  (6  +  1/4) 

2 

and  the  lower  line  is  a  -  8,  where  a  -  (fl  (y  k  C)  is  normalized 
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strength  of  the  wiggler  field  and  «  -  u?/{2y*  B2  k2C2)  is  normalized 
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beam  plasma  oscillations. 
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A 

The  solid,  dashed-dotted  and  dotted  lines  shows  as  a  function  of 
ax  for  or  »  10,  20  and  30,  respecively. 


-20  -10  0  10  20  30  40  50  60  70 

(a-Aa>)r 


Fig.  5  Graph  of  normalized  gain  G  u  a  function  of  a_x  for  <*t  -  20, 
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Appendix  A 

We  are  interested  In  relating  the  parameters  p+,  p_,  9+  and  8_  in  Eq. 

(56)  with  the  initial  conditions  of  the  electron  beam  Injected  into  quadrupole 
periodic  magnetic  field.  Electron  position  is  described  in  cylindrical 
coordinates  by  its  radius  R  and  angle  8,  radial  velocity  Vr  and  angular 
velocity  Vg.  From  Eq.  (56)  follows  that 


x^°\o)  »  p+cos8+  +  p_cos8_  ■  R  co 8  8 

y^(0)  ■  p+sin8+  -  p_sin0_  *  R  sin  0 

v£o)(0)  -  -a(p+sin8+  +  p_sin0_)  -  VRcos0  -  V0sin0 

Vy°^(0)  -  a( p+cos8+  -  p_cos8_)  »  VRsin8  +  V0cos8 

then 

V  V 

p  CO80  T  [(R  +  -^)  cos 8  +  -£•  sin8  ] 
r  t  ■  4  a  a 

V  V 

p  cos8  -  -7  l(R - -)  cos0 - —  sin0} 

■  ~  *  a  a 

1  VB  VD 

p  sin8  ■  M(R  +  — )  sin8 - -  cos8] 

+  +  2  a  a 

V  V 

p  sin8  -  -  -~  [  (R - -)  sin8  +  —  cos8  ] . 

~2a  a 


(A— 1 ) 


(A-2) 


Introduce  a  following  notation 
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+  —  *  u,  cosy, 
ail 


a - ^sin^ 


(A-3) 


R  -  —  *  U2cos^2 


—  -  U  sirnj>0 
a  2  2 


(A-4) 


from  (A-2)  we  obtain 


p+CO80+  ■  jDj  C08  (0  +  ij^) 


p+sin0+  -  y  Ux  sin  (0  +  i^) 


p_cos0_  -  y  U2  cos  +  ^2^ 


p_sin0_  “  “  7  U2  sin  +  V 2^ 


(A-5) 


with  constraint 


U2  sin  ^  "  U1  8in  V 


(A-6) 


we  can  express  P+,  P_,  0+  and  0_  in  terms  of  R»  0,  and  ♦  j,  namely  from 


(A-5)  follows 


